Radio frequency association of heteronuclear Feshbach molecules 
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We present a detailed analysis of the production efficiency of weakly bound heteronuclear ^"K *^Rb 
Feshbach molecules using radio frequency association in a harmonic trap. The efficiency was mea- 
sured in a wide range of temperatures, binding energies and radio frequencies. A comprehensive 
analytical model is presented, explaining the observed asymmetric spectra and achieving good quan- 
titative agreement with the measured production rates. This model provides a deep understanding 
of the molecule association process and paves the way for future experiments which rely on Feshbach 
molecules e.g. for the production of deeply bound molecules. 



The creation and investigation of ultracold, trapped, 
deeply bound molecules has been a major research goal 
for a number of years. In particular polar molecules with 
long range dipolar interaction [ij offer a number of fas- 
cinating experimental avenues such as the investigation 
of quantum phases 0, methods for quantum computa- 
tion [3||or precision measurements of fundamental con- 
stants [J|. Recently, a number of approaches have made 
significant progress towards that goal. Polar molecules 
have been slowed down and trapped using Stark decel- 
eration in spatially modulated electric fields 0, @, Q- 
Other approaches are provided by cooling ground state 
molecules with a buffer gas [§| and by photoassocia- 
tion in magneto-optical traps, where the lowest vibra- 
tional level in ultracold K2 H, KRb [l^, RbCs [lH, 
Cs2 [3, and LiCs has been reached. However, 
these approaches only offer limited state control or tem- 
peratures in the mK regime. On the other hand, the 
understanding and control of interactions in ultracold 
gases has made dramatic progress in the past few years. 
The cooling procedure for two atomic species as well as 
the heteronuclear interaction potentials have been stud- 
ied in detail 14, [3, [Til . [itI, e speciall jy in the case of 
"OR- and *^'^Rb-mixtures [U lH, • This develop- 

ment now allows for the association of weakly bound het- 
eronuclear molecules in well defined quantum states using 
magnetically tunable Feshbach resonances [S^, 23, 24 1. 
A detailed understanding of the association process of 
these molecules is of central importance. These Feshbach 
molecules provide an ideal starti ng p oint for a further 
transfer into deeply bound states [25, 2^, 27 1. Recently, 
the de-excitation of these molecules to the rovibrational 
ground state has also been demonstrated [2^. 

So far, the association of homonuclear molecules has 
been described with a semi-classical numerical model 
[2^, which does not predict absolute conversion rates. 
The model assumes adiabatic conversion of selected atom 
pairs to molecules. It has also been applied to het 
eronuclear mixtures 12 



24| . However, in this case the 



magnetic fields for the production of molecules 3^, 31 1 . 



approach is only valid if the conversion speed is small 
enough to assure adiabaticity, while it has to be large 
enough such that the atom pairs do not drift apart. Fur- 
ther theoretical studies focus on the use of modulated 
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FIG. 1: Number of associated molecules versus radio fre- 
quency for various binding energies. The solid lines show the 
theoretical model. All measurements were taken at a temper- 
ature of 730 nK. The duration of the radio frequency pulses 
is indicated in the legend. The inset illustrates the association 
process including the relevant energies and the wave functions 
of molecule and atom pair. 

In this paper, we present a complementary analytical 
model for radio frequency association which yields abso- 
lute molecular production rates. The model is compared 
to radio frequency association spectra obtained in an ul- 
tracold mixture of '^^Rb and "^^K atoms. It explains 
the measured asymmetric molecular association spectra 
depending on the energy distribution of the atoms in 
the trap and thus provides a deep understanding of the 
molecule association process. 

In our approach, the number of molecules is calculated 
by time dependent perturbation theory assuming a Gaus- 
sian time dependence of the coupling term. The treat- 
ment is analogous to the derivation of Fermi's Golden 
Rule, which has been successf ully used to describe ra- 
dio frequency dissociation 32, l33[ . In the case of radio 
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frequency association, the number of colliding pairs per 
energy interval h{er) has to be included, where is the 
energy of the relative motion. In the limit of small deple- 
tion, the perturbative approach for short radio frequency 
pulses with a Gaussian shaped amplitude of duration r 
yields the number of molecules. 
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Here, Erf = hi> is the energy corresponding to the 
radio frequency and H is the Rabi frequency of the 
atomic transition with energy Eq (sec inset of Fig. [1]). 
Ef, denotes the molecular binding energy and 



Ff{er) 



(2) 



is the Franck-Condon factor between the wave function 
'J'j^ (r) of the colliding atom pair and the boimd molecu- 
lar wave function $,„(r). 

To evaluate equation ([1]), we start with the occupa- 
tion density h{er) by considering the distributions of the 
two atomic species in the harmonic trap. The Hamilto- 
nian of two noninteracting atoms a and b describes two 
uncoupled three-dimensional harmonic oscillators 



number of possible atom pairs per state and total energy 
et = ea + Eb is simply the product of the two single occu- 
pation densities /p(eO = NaNi,i^)^exp{-et/ikBT)), 

where uj = (Oi'^i)^^'^- Since the two Hamiltoni- 
ans Hi and H2 describe the same system, degenerate 
states in the two-atoms-representation have to map the 
corresponding degenerate eigenstates in the atom-pair- 
representation homogeneously. Thus, the number of 
pairs per state /p(et) has to be equal for the two rep- 
resentations. In order to obtain the total number per 
energy, the density of states has to be included, which is 
the product of the density of states of the relative and 
the collective motion. The three-dimensional harmonic 
oscillator of the center-of-mass motion has a density of 
states gcm{£cm) = 2(hEy^ • O'^ly atom pairs colliding with 
vanishing angular momentum contribute to the forma- 
tion of s-wave Feshbach molecules. Therefore, the den- 
sity of s-wave states per relative kinetic energy is given 
bygr(er) = 2fe. 

For the molecule production, only the energy of the 
relative motion is relevant, since the molecule follows 
the same ccntcr-of-mass motion as the atom pair. The 
number of molecules per energy in the relative motion 
can be calculated by integrating over all possible ccntcr- 
of-mass energy states. 
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where Pa/b denotes the momenta of the atoms, Xa/b,i 
the ith component of their positions, m^/b their masses 
and ^0Ja/b,i the ith component of their trapping fre- 
quencies. One can introduce center-of-mass coordinates 
as usually done in free space, where r,p,^ describe po- 
sition, momentum and mass of the relative motion and 
i?, P and M the motion of the center of mass [3l| . 

The Hamiltonian then reads 
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where uji — {ujaA'-^b,iY^'^ is the mean trapping fre- 
quency. This Hamiltonian again contains two harmonic 
oscillators, one for the relative and one for the collective 
motion. The residual term He couples relative and col- 
lective motion and is proportional to (wa,i — ^b,i)- If the 
trapping frequencies are similar, this coupling is small 
and will be neglected for the molecule production on 
short timescales. 

The number of atoms per state and energy e in 
a single-atom harmonic oscillator is assumed to be 
Maxwell-Boltzmann distributed in the thermodynamic 



limit /(e) = N 



exp{—e/{kBT)). Therefore, the 
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For the complete evaluation of eq. ([T]), we now evalu- 
ate the Franck-Condon factor Ff{er). In the case of small 
interatomic distance, the wave function of the atom pair 
matches the wave function of a free atom pair ex- 
cept for a factor of huj for correct normalization. To 
calculate the Franck-Condon factor, it is crucial to in- 
clude the scattering length of the colliding atom pair a', 
since the molecule production probes the innermost part 
of the scattering wave function. Nevertheless, it is valid 
to neglect the influence of the scattering length on the 
energy distribution at the relevant densities. Hence, we 
adopt the Franck-Condon factor from Ref. [ssj . 
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where Eb is the binding energy of the formed molecule 
and according to Ref. [33| , a and i?^ are defined implicitly 
by 
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The factor xi^b) reflects that the molecular wave func- 
tion consists of two parts, the closed channel and the open 
channel fraction. Only the latter contributes significantly 
to the Franck-Condon factor. It can be approximated [2J| 
by x{Eb) = 1 - dEt/dB, where B denotes the 

magnetic field and A/i is the difference in magnetic mo- ^ 
ment of the molecular and the atomic state. In the vicin- 3 
ity of the Feshbach resonance (sign(abg) b-Bq ^ 
can be approximated using ([7]) and the scattering length uf 



a{B)=abg{l-AB/{B~Bo)). 



(8) 



Here, abg denotes the background scattering length, Bq 
and AB are the position and the width of the resonance. 
In total, one obtains 



xiEb) = 1 - 



AfiAB ji abg 



(1 4- kabg 
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where fc^ = Eb/2^. 

Now, all parts of equation ([T|) are derived. In the fol- 
lowing, we directly fit the theoretical prediction to our 
experimental data, allowing only two free parameters: 
The binding energy Eb and a scale factor A to quantify 
deviations from our model. For the atom pair scatter- 
ing length, the values abg = 9.88 nm, a' = 9.10 nm and 
= 2.32/iR are adopted from a coupled channel cal- 
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We used the following experimental sequence to pro- 
duce heteronuclear ""^K-^^Rb-molecules close to a Fesh- 
bach resonance located at 547 G and to measure their 
production efficiency. First, a mixture of 2 x 10^ ^°K- 
Atoms and 1 x 10^ ^^Rb-atoms at a temperature of 2 fiK 
is prepared in a crossed beam optical dipole trap at a 
wavelength of 1064 nm. We employ rapid adiabatic pas- 
sages to transfer the ^^Rb atoms to the |1, 1) state and 
the atoms to the 1 9/2, —7/2) state at homogeneous 
magnetic fields of 10 G and 20 G. Then the magnetic field 
is increased to values below the Feshbach resonance. The 
sum of magnetic field noise, field inhomogeneity and shot- 
to-shot fluctuations is less than 16 mG, accomplished by 
a mechanical transport to the geometric center of the 
Feshbach coils [s^ . The power of the dipole trap is low- 
ered to mean trapping frequencies of 244 Hz (335 Hz) for 
^^Rb (*°K), yielding samples of 5 x 10'^ '^°K atoms and 
2.5 x 10'^ ^''Rb atoms at a temperature of 730 riK. Fi- 
nally, we apply a Gaussian shaped radio frequency pulse 
with a width of t = 12.5 to 31.25 fis and a peak atomic 
Rabi frequency of f2 = 27r 45 kHz to produce up to 5x 10"* 
weakly bound molecules (see inset of Fig. [1]). The dura- 
tion of the radio frequency pulse is chosen short enough 
to ensure that the production efficiency is still in the lin- 
ear regime, i.e. the loss of molecules due to interaction 
with residual atoms can be neglected [s^ . 

Immediately after the molecule association, the trap- 
ping potential is switched off and the number of associ- 
ated molecules is measured directly using resonant light 



for atoms in the state |9/2, —9/2) after 2 ms of ballis- 
tic expansion. In two subsequent absorption images, we 
measure the remaining number of ""^"K and ®^Rb atoms 
in the states |9/2,-7/2) and |1,1). 
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FIG. 2: The molecular binding energies as a function of 
magnetic field are shown in the top frame. A value of 
Bo = 546.618(6) G is obtained for the position of the reso- 
nance (dotted line). The bottom shows the obtained number 
of molecules per pulse length. 

In a first set of experiments, we have studied the effi- 
ciency of molecule production for six magnetic fields be- 
tween 545.73 G and 546.19 G. For each magnetic field, 
the atomic transition frequency of the *'^K atoms was 
measured as a reference for the association spectra. Fig- 
ure [T] shows the absolute number of KRb molecules ver- 
sus radio frequency detuning for various magnetic fields. 
The asymmetry of the association spectra is clearly vis- 
ible and reflects the relative kinetic energy distribution 
of the associated atom pairs. The shape of the spectra 
is very well reproduced by the model described above, 
where only the binding energy Eb and the scaling factor 
A were used as fit parameters. 

In Figure [2l the extracted binding energies and 
molecule production rates are plotted as a function of 
the magnetic field strength. A fit of eq. ([7]) to the data 
yields Bq = 546.618(6) G (calibration uncertainty 5 mG) 
for the position and AB = 3.04(2) G for the width of the 
resonance. This combination of experimental measure- 
ment and theoretical analysis represents the most precise 
characterization of a heteronuclear resonance in an opti- 
cal dipole trap, consistent with Ref. [l^l- The width of 
the resonance enters the model through eq. Q and a 
self-consistent description was obtained by iteration. 

Figure |4] (a) shows the results for A which represents 
the total deviation from the theoretical model. It shows 
that this model quantitatively predicts the heteronuclear 
molecule yield with good accuracy. The comparably 
small overall deviation by a factor 2 can be explained as 
follows: Primarily, coUisional relaxation leads to a loss 
of molecules during the expansion time [s^ . In addition 
to the overall deviation, the molecule production rate 
drops for small binding energies close to the Feshbach 
resonance, since the coUisional relaxation is amplified by 
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the resonantly enhanced scattering length. Moreover, the 
model depends on the temperature, the trapping frequen- 
cies and the product of ^^Rb and "^^K atom numbers, 
thus slight miscalibrations strongly affect the results. 
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FIG. 3: Number of associated molecules versus radio fre- 
quency for various temperatures. The solid lines show the 
theoretical model. All measurements were taken at a fixed 
magnetic field of 545.994 G. 

In a second set of experiments, we have investigated 
the dependence of the molecule production on the tem- 
perature of the atomic mixture. By lowering the intensity 
of the two trapping beams, ^°K-®^Rb-mixtures in a tem- 
perature range from 1.1 fj,K down to 250 nK with cor- 
responding mean trapping frequencies of 240 to 170 Hz 
(340 to 230 Hz) for 87Rb (40k) were produced. The 
association spectra are shown in Fig. [3] with theoreti- 
cal fits to the data, again using only binding energy Eb 
and the scaling factor A as free parameters. Once more, 
the theoretical distribution functions match the exper- 
imental findings very well. Since the magnetic field is 
kept constant at 545.994 G, the resulting binding ener- 
gies of all independent fits yield a value of 127.6(6) kHz, 
as expected. Figure [J] (b) shows the scaling factor A as 
a function of temperature. It drops slowly with decreas- 
ing temperature, since the different quantum statistics 
of the two species start to play a growing role. While 
the bosons accumulate in the lowest energy levels, the 
fermions obey the Pauli principle. This yields a smaller 
overlap between the clouds and thus a smaller molecule 
production efficiency. Additionally, the internuclear in- 
teraction has been neglected in the occupation densities 
of the atomic clouds. This repulsive interaction leads to 
smaller effective trapping frequencies and accordingly to 
a smaller molecule yield. 

In conclusion, we have measured the production ef- 
ficiency of heteronuclear Feshbach molecules in a wide 
parameter range of temperatures, binding energies and 
radio frequencies. A comprehensive analytical model is 
shown to provide a detailed understanding of the ob- 
served spectra and a good quantitative estimate of the 
absolute production rates. It is of particular relevance 
for ongoing experiments with heteronuclear Fermi-Fermi 
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FIG. 4: Scaling factor versus binding energy (a) and temper- 
ature (b). The solid lines are guides the eye. 



mixtures [l^, 17 1. This work thus paves the way for 
further experiments using heteronuclear molecules as a 
starting point for the production of ultracold polar sam- 
ples. 

It will be especially interesting to include the influ- 
ence of quantum statistics in the model. A numerical 
approach will be subject of future work. 

We acknowledge support from the Centre for Quantum 
Engineering and Space-Time Research QUEST and from 
the Deutsche Forschungsgemeinschaft (SFB 407, Euro- 
pean Graduate College Interference and Quantum Ap- 
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Note added: Recently, we have become aware of a com- 
plementary work with a bosonic heteronuclear mixture of 
^^K and M. 
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